Measurement of exclusive D meson decays to 77 and 77' final states and SU(3) 

amplitude analysis 
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Using 281 pb _1 of data collected with the CLEO-c detector, we present new measurements of 
Cabibbo-suppressed decays of D° and D + mesons to n and n' final states. We make first observations 
of D° — > r/'ir , nn, r\rj , and t]-k + -k~ , and find evidence for D + — » 7?7r + 7r , D + — > ti'ty + tt° and D° — > 
n'-K + Tr~ . We also report on improved measurements of D° — > niv , D + — > rj7v + and D + — > n'-R + . 
Using the measured two-body Cabibbo-suppressed decays, we extract amplitudes for specific flavor 
topologies and compare them to those from Cabibbo-favored decays. 

PACS numbers: 13.25.Ft 



Charm decays provide a laboratory for the study of the weak and strong interactions. Because of their sim- 
plicity, two-body hadronic decays provide experimental input for understanding strong final state interaction 
effects in heavy meson decays [J [H, H, EL H, HI- By studying these decays, one can determine the magnitude 
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and relative phases of contributing isospin or topological decay amplitudes. To describe these decays, partic- 
ularly for D and B mesons, one often invokes symmetries, such as SU(3) flavor symmetry, which assumes the 
masses of the u, d and s quarks are equal (and implicitly assumed to be small). Small differences in decay 
rates are then attributed to SU(3) symmetry-breaking effects. Typically, such symmetry-breaking effects are 
at the level of 20% for D meson decays, and are expected to decrease as the parent meson mass increases, 
e.g., B decays. In the SU(3) flavor quark diagrammatic approach [J], two-body decays are decomposed into 
contributions from tree (T), color-suppressed (C), annihilation (A) and exchange (E) diagrams. Penguin 
diagrams are expected to be negligible in charm decays and are not discussed any further. Additional dis- 
connected diagrams, such as the singlet-exchange (SE) and singlet-annihilation (SA) may also contribute to 
final states that include r\ and rf via their coupling to the SU(3) singlet component of these mesons. These 
quark diagrams are shown in Fig. [TJ 




(c) Exchange (E). (d) Annihilation (A). 
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(e) Singlet Exchange (SE). 



FIG. 1: Feynman diagrams for various flavor topologies used to describe D — > PP decays. Here, we use Cabibbo- 
favored diagrams as an example, for T, C, E, and Cabibbo-suppressed decays for the A and SE diagrams. 



An analysis of D — > PP Cabibbo-favored (CF) decays has been carried out [7j], and it was shown that 
their branching fractions are well-described within the SU(3) flavor topology approach. In this report, we 
measure branching fractions of singly Cabibbo-suppressed (CS) decays and analyze them within the same 
framework Q. 

This analysis utilizes 281 pb _1 of data collected on the -0(3770) resonance at the Cornell Electron Storage 
Ring. This energy is just above threshold for production of DD pairs, and thus no additional particles 
accompany the pair. The decay products of the D mesons are reconstructed using the CLEO-c detector, 
which is a general purpose solenoidal detector. The detector includes a tracking system for measuring 
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momenta and specific ionization (dE/dx) of charged particles, a Ring Imaging Cherenkov Counter (RICH) 
for particle identification, and a Csl calorimeter (CC) for detection of electromagnetic showers. The CLEO-c 
detector is described in detail elsewhere @. 

We reconstruct D meson candidates in the following singly Cabibbo-suppressed decay modes: D + — > r/7r + , 
rfir + , r/Tr + ir°, i]'ir + ir a , and D° — > i]ir° , ij'ir , r/rj, rjrj' , ?/7r + 7r _ and r)'ir + ir~. Unless otherwise noted, charge 
conjugate final states are assumed throughout. 

Charged particles are reconstructed using the tracking system and are required to pass a set of standard 
selection criteria [l(J. Charged pions are identified using a \ 2 discriminant based on dE/dx and RICH 
information. First, we require that the measured dE/dx is within three standard deviations (a) of the 
expected value for a pion at the given momentum. We then define: 



A x 2 = x 2 k-xI +LL K -LK, (1) 

where Xk(-k) is the difference in the measured and expected dE/dx for the kaon (pion) hypothesis, normalized 
by its uncertainty, and LL^m is the negative log-likelihood for the kaon (pion) hypothesis obtained from 
the RICH. The RICH information is used only for charged particles with momentum above 0.7 GeV/c and 
possessing at least three associated Cherenkov photons. We require A\ 2 > 0, which results in an efficiency 
of about 95% and a fake rate of no more than 2%. We construct ir° — ► 77 and 77 — > 77 candidates using 
reconstructed showers in the calorimeter that have energy E 1 > 30 MeV [E 1 > 50 MeV for 77 — > 77) and no 
charged track that is projected within the vicinity of the shower. For each candidate, we require the mass pull, 
Cm = (M lec — Mp) I a M rcc < 3, where M rec is the reconstructed invariant mass, Mp is the parent particle mass 
(here ir° or 77), and <TM xec is the uncertainty in the invariant mass of the candidate («5-7 MeV/c 2 ). Candidate 
77 — > 7T + 7r~7r° decays are also formed using the previously discussed selection criteria on charged and neutral 
pions. To improve the resolution on the tt° or 77 momentum, the decays are kinematically constrained to 
the known ir° and 77 meson masses. The improved four-momenta are then used in subsequent analysis. We 
select 77' — > 777r + 7r~ candidates by requiring the mass difference, 402 < — M n < 418 MeV/c 2 . 

For D decay modes with only a single 77, there is minimal benefit to including the 77 
because of the lower branching fraction and detection efficiency. We therefore only use 77 — * 77. However, for 
D — > 7777 and D° — ► 7777', there is a significant gain in statistical power by allowing one 77 — ► 7r + 7r~7r° decay 
in the final state. To reduce combinatoric background in three-body decays we veto 77 — > 77 candidates 
that share a photon with any tt° with \<jm\ < 3.0. To ameliorate backgrounds from Cabibbo-favored modes 
containing a K s we require oppositely charged pion pairs in the final state to be outside the invariant mass 
interval from 475-520 MeV/c 2 . 

Signal D candidates are formed and required to have an energy, Ed, consistent with the beam energy, 
-E-beam, by requiring AE = .Ed - -Eticam is consistent with zero. The mode-by-mode AE selection requirements 
are determined from Monte Carlo simulation, and are shown in Table U Efficiencies for these decays after 
all analysis requirements are determined from Monte Carlo simulation of these decays [l5|, [l6|, [TtJ and are 
also shown in Table U 

For candidates passing these selection criteria, we compute the beam-constrained mass, M^ c — 
V^beam — IP-Dp' where po is the momentum of the D candidate. Substituting -Ebcam for the candidate 
energy improves the mass resolution by about a factor of two. Figure [2] shows the M^ c distributions for 
the two-body singly Cabibbo-suppressed decays (a) D + — > 7777+, (b) D° — > i]ir°, (c) D + — > 77'^+, and 
(d) D° — > 77'7r°, (e) D° — > 7777, and (f) D° — * 77'?7. Prominent peaks are observed for all six decay modes. 
The points show the signal candidates in data and the curves are fits based on maximizing the likelihood 
of the probability density function, which is given by the sum of an ARGUS threshold function [l2l ] and an 
asymmetric Gaussian signal shape (cbal) [l3j. The ARGUS shape parameters are extracted by fitting Mb c 
distributions obtained from the AE sideband regions in data. The high-mass tail in the M^ c spectrum 
results from initial state radiation (ISR) and is modeled by a power-law tail in the CBAL line-shape. The 
signal shape parameters (mean, width and tail parameters) are determined from, and fixed to the values 
obtained from Monte Carlo simulation. We obtain yields of 1033±42 for D+ -> 7771-+, 160±24 for D° -> 7777°, 
352 ± 20 for D+ -> 77^+, 50 ± 9 for D° -> t/tt , 255 ± 22 for D° -> 7777, and 46 ± 9 events for D° -> r(r\. 
These are first observations of D° — > ?7'7r , D° — > 7777 and D° — > 77' 77, and correspond to signal significances of 
7.2, 14.8. and 6.7, respectively. The significances arc obtained from the differences in log-likelihood values 
with and without (signal yield set to zero) the signal component i.e. y2(A logL) . The yields are shown 
in Table HU For D —> r\r\ and D° —> rjrj' , we also show yields when the two 77 decay channels are analyzed 
independently; both are consistent with fits to the sum of both decay channels. 
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TABLE I: Summary of mode-dependent AE selection requirements for modes under consideration. Also shown are 
the efficiencies from signal Monte Carlo simulation. For D° — > 7777 and D° — > rjrj' , we indicate the decay modes of the 
77's in parentheses. The quoted uncertainties are statistical and systematic, respectively. 



Mode 




AE Range (MeV) 


Efficiency (%) 


D+ - 


7771" + 


[-28, 25] 


46.7 ±0.5 ±1.7 


D+ - 




[-19, 18] 


27.7 ±0.5 ± 1.1 


D° -» 


777T 


[-45, 34] 


30.6 ±0.5 ±1.3 


D° -» 




[-38, 32] 


17.1 ±0.4 ±0.7 


D° -» 


w (77) (77) 


[-33, 30] 


28.9 ±0.4 ±2.2 


D° -» 


W (77)(7r + 7r _ 7r°) 


[-26, 25] 


16.5 ±0.4 ±0.7 


23° -» 


W (77) (77) 


[-27, 23] 


15.9 ±0.4 ±1.2 


D° 


7777' (77)(7T + 7r~7r°) 


[-23, 20] 


8.4 ±0.3 ±0.4 


Z?° -» 




[-22, 19] 


29.1 ±0.5 ±1.1 


D+ - 


r]n + 7r° 


[-29, 24] 


16.7 ±0.4 ±0.7 


D° -» 


Tj'n + n~ 


[-19, 17] 


13.3 ±0.3 ±0.5 


D+ - 




[-22, 20] 


7.3 ± 0.3 ± 0.3 



Peaking backgrounds from non-resonant final states, such as D° — > r)ir + ir~ ir° , D° — > r)rjir + ir~ , etc, which 
may contribute to the D° — > r\r\ and D° —> r]r)' decays, respectively, are highly suppressed because the 
7r + 7T~7r° and r]Tr + ir~ invariant masses must be consistent with the 77 and 77' masses, respectively. Moreover, 
these backgrounds are also Cabibbo suppressed, and thus we expect them to be very small or negligible. 
This supposition is checked by selecting candidates from the sideband regions, 7 < \<tm\ < 10 and 8 < 
M lt + V - V — M n — 410 < 16 McV/c 2 , and repeating the analysis. No evidence of peaking backgrounds are 
found. 

Figure [3] shows the M\, c distributions for the three-body Cabibbo-suppressed decays, (a) D Q — > r)iT + Tr~, 
(b) D + — > 777r + 7r°, (c) D° — > ri'ir + ir~, and (d) D + — > rj'Tt + ir a . We obtain signal yields (significances) of 
258 ±32 (9.0) for D° -> rjn+TT-, 147 ±34 (4.5) for D+ -> tjtt+tt , 21 ±8 (3.2) for D° -> rj'n+n^, and 33 ±9 
(4.2) for D + — > 77'7r + 7r°, respectively. We thus establish the Z?° — * ?77r + 7r _ decay, and provide first evidence 
for the other three-body decay modes. 

The branching fractions are computed from: 



where iV S i g is the number of signal events, N d q is the number of DD pairs produced, and the acceptance 
A = 5^e<fJ£. Here, the sum is over the product of efficiency (ej) and rf 1 ' submode branching fractions, 
for the measurements indicated in Table HT1 The numbers, N D o£,o = (1.031 ± 0.015) x 10 6 and N D + D - = 
(0.819 ± 0.013) x 10 6 , are determined from an independent measurement of Cabibbo-favored D hadronic 
branching fractions [H| . The yields and branching fractions are tabulated in Table [TXJ For D° — » 7777 and 
D° — * 7777' we also show that the two decay modes yield consistent branching fractions. 

The inclusive rates for D — > -q^X have been recently measured [14], and the exclusive D° — > tjX, 
D° — > r]'X and D + — * 77X modes measured here comprise about 10% of their respective total inclusive rates. 
In contrast, B(D + — > t]'tt + ) ± B(D + — > ?7'7r + 7r ) ~ 0.6% accounts for about 60% of the total inclusive rate. 

A number of systematic uncertainties have been considered. Efficiencies, as determined from Monte Carlo 
simulation, are subject to uncertainties due to finite statistics (1-2%), modeling of the underlying physics, 
and modeling of the detector response. The underlying physics that induces uncertainty in the efficiencies 
includes modeling of final state radiation (1%) and resonant substructure. The latter is applicable only to the 
three-body decay modes and is determined by comparing efficiencies determined using a phase-space decay 
with those obtained using intermediate resonances, such as Tj^p, ao(980)7r, or ao(1450)7r. We take the largest 
fractional difference in efficiency as the associated uncertainty. The values range from 3% for D° — ► t]it + -k^ to 
12% for D + — > i]'ir + Tr a . Charged and neutral particle reconstruction and identification has been extensively 
studied using a missing-mass technique [lOj], and we find that the efficiencies in data are consistent with, 
or slightly lower than simulated efficiencies. We thus correct the simulated efficiencies as follows. The pion 
tracking efficiency and particle identification corrections are 1.000 ± 0.003 and 0.9950 ± 0.0025, respectively. 



5 



TABLE II: Summary of yields and branching fraction measurements, as discussed in the text. For D — > rjrj and 
D — > rjTj' , we also show the individual results obtained from the two r\ submodes. The first uncertainty is statistical 
and the second is systematic. Where measurements are available, results are compared to the PDG. 



Mode 


Yield 


Branching Fraction 


PDG[11] 

(10- 4 ) 








(10- 4 ) 


D + 


— > 717Y^~ 


1033 ± 42 


34.3 ±1.4 ±1.7 


35.0 ±3.2 


D+ 


/ + 
—* 7] 7T^ 


352 ± 20 


44.2 ±2.5 ±2.9 


53 ± 11 


D° 




156 ± 24 


6.4 ± 1.0 ±0.4 


5.6 ± 1.4 


D° 


l 

— > r? 7r 


50 ±9 


8.1 ±1.5 ±0.6 




D° 




255 ± 22 


16.7 ± 1.4 ± 1.3 






(77) (77) 


141 ± 17 


15.3±1.8(stat.) 




D° 


(77)( 7r+7r_7r °) 


115 ± 13 


19.0 ± 2.2(stat.) 




-» 777/ 


46 ±9 


12.6 ±2.5 ±1.1 






(77) (77) 


33 ±8 


14.8 ± 3.3(stat.) 






(77)( 7r+7r_7r °) 


14 ±5 


10.5 ± 3.5(stat.) 




D° 


— > rjn + n~ 


257 ± 32 


10.9 ±1.3 ±0.9 


< 19 


D+ 


— > r)TT + TT° 


149 ± 34 


13.8 ±3.1 ±1.6 




D° 


— > rj'n + n~ 


21 ±8 


4.5 ±1.6 ±0.5 




D+ 


t/tt+tt 


33 ±9 


15.7 ±4.3 ±2.5 





For each tt° — ► 77 (77 — * 77), we adjust the efficiencies by 0.961 ± 0.020 (0.943 ± 0.035). These corrections 
and uncertainties are included in the efficiencies shown in Table [U 

Additional sources of uncertainty arise from the candidate selection requirements, namely the Kg veto, 
rl mass window, r\ — > 77 veto (for photons also used in a 7r°), and the restricted AE range. We estimate 
the uncertainty introduced from these requirements by increasing the window size and taking the fractional 
difference in the efficiency-corrected yield between the nominal and the larger window size. For the K s 
rejection, we increase the veto region by ±50%, from which we find an uncertainty of 2%. For the 77' selection, 
we broaden the selection window to cover the range from 395-423 MeV/c 2 , and find an associated uncertainty 
of 1%. The uncertainty introduced from the 77 — > 77 veto is only applicable to three-body decay modes. It 
is estimated by comparing (data and MC simulation) the fraction of signal events that pass the veto with 
respect to no veto in the two-body decay modes. The ratio, averaged over several modes, is 1.01 ± 0.02, and 
we conservatively assign a 5% systematic error to this source. The AE systematic uncertainty is obtained by 
increasing the AE signal window by one unit of the Gaussian width. The corresponding uncertainty ranges 
from 2% for D + — ► r]ir + to 4% for D° — ► r]ir . Uncertainty in the signal yield receives contributions from both 
the signal and background shape parameterizations. Both are studied by varying the shape parameters one 
at a time by ± 1 standard deviation and adding the resulting changes in yield in quadrature. The associated 
uncertainties range from 1% (for D — > rjir + ) to 4% (for D° — > r/7r + 7r _ ) for the signal shape uncertainty and 
0.5% (for D — > t]'tt + ) to 4.0% (for D° — > rfir + ir~) for the background shape and normalization. Uncertainty 
due to multiple candidates in an event is quantified by computing the average number of candidates per 
event within ±5 MeV of the known D mass for both simulation and data. We find good agreement between 
data and simulation, except for final states with low momentum 7r°'s, where the data has a slightly larger 
rate of multiple candidates. We take half the difference between the average number of multiple candidates 
between data and simulation as a correction to the branching fraction, and assign 100% uncertainty to it, 
leading to the following corrections: 1.013 ± 0.013 for D + — > 7?7r + 7r°, 1.015 ± 0.015 for D + — > r/7r + 7r , and 
1.03 ± 0.03 for D° — > (77 — > 77X77 — > ir + ir~ir°). The rest of the modes are consistent with unity at the level 
of 1%, which we assign as a systematic uncertainty. The branching fractions for 77 — > 77, 77 — > Tr + ir~ir° 
and 77' — * ?77r + 7T~ are uncertain by 0.6%, 1.5% and 3.1%, respectively, for each such decay in the final state. 
Lastly, the number of D°D° (D + D~) events in our data sample has an uncertainty of 1.5% (1.6%) [10(. The 
systematic uncertainties for all modes under consideration are summarized in Tables IIIII and are included in 
the branching fraction measurements in Table HT1 

We have searched for intermediate resonances in D° — > i]ir + ir~ . Figure 2] shows the sideband-subtracted 
777T + and 7r + 7r~ invariant mass distributions. Surprisingly, there are no significant contributions from either 
•qp or ao(980)?T + . Overlayed on the data (points) is a Monte Carlo simulation where a phase space model 
is used. We find that decay is well-modeled by three-body phase space. 
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1.83 1.86 1.89 1.83 1.86 1.89 

M bc (GeV/c*) 

FIG. 2: Distribution of Mb c for the two-body Cabibbo-suppressed decay modes: (a) D + — > r?7r + , (b) D + — > r]'ir + , (c) 
D — > t/tt , (d) D° — ► ry'vr , (e) _D° — > 7/77, and (f) Z?° — > 7777'. The superimposed curve is a fit to the data as described 
in the text. 

We now apply an SU(3) diagrammatic analysis to test the validity of the SU(3) flavor decomposition 
approach to charm meson decays 0, 0, i]. Two-body D meson decays can be described using an SU(3) 
diagrammatic approach in terms of an external tree diagram, T; a color-suppressed tree diagram, C; an 
exchange diagram, E; an annihilation diagram, A; and a singlet exchange diagram, SE. The SE contribution 
represents the matrix element that produces an 77 or 7/ through its coupling to the SU(3) singlet portion of 
these mesons. Such contributions are OZI-suppressed, and thus expected to be small. These five diagrams 
are shown in Fig. [TJ 

We first update the fitted diagrammatic amplitudes obtained from Cabibbo- favored decays Q . Table IIVI 
shows the SU(3) representations, the measured branching fractions, decay momentum (p*), and invariant 
amplitudes for Cabibbo-favored D — * PP decays. The lifetime of the D mesons and the D — > K°(r),i]') 
branching fractions are taken from the PDG whereas the other six branching fractions are taken from 
recent CLEO measurements EE Gjl • The eight branching fraction measurements are fit to four topological 
amplitudes and three relative strong phases. The phases of C, E and A are measured with respect to T, 
which is taken to be real. The phase, 5ab is the angle subtended from amplitude B to amplitude A. The 
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1.83 1.86 1.89 1.83 1.86 1.89 



M bc (GeV/c*) 

FIG. 3: Distribution of M^ c for the three-body Cabibbo-suppressed decay modes: (a) D° — » t]ty + it~, (b) D + — > 
t)-k + -k , (c) D° — > ?7'7r + 7r _ , and (d) D + — > r;'7r + 7r°. The superimposed curve is a fit to the data as described in the 
text. 

fitted amplitudes (in units of 1CP 6 GeV) and relative strong phases are found to be: 



T = (2.78 ±0.13), (3a) 

C = (2.04±0.17)e ,( " 151±2) °, (3b) 

E = (1.68±0.12)e l(117±4) °, (3c) 

A = (0.54±0.37)e l( - 64± « 2) °. (3d) 



The fit x 2 is 0.65 for 1 degree of freedom, indicating that the OZI-suppressed SE and SA contributions are 
not needed to describe the branching fraction measurements in Cabibbo-favored decays. These results are 
consistent with those obtained in Ref. [|| and the fitted branching fractions (B&t in Table Hv|) are in good 
agreement with the experimental values. 

We now analyze the Cabibbo-suppressed decays using the same framework. Using our measured branching 
fractions for D° — ► 7771"°, D° — > rfir , D° — » r\r\ and D° — > r\rj , we extract C and E', the analog of C and E 
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TABLE III: Systematic uncertainties for signal D decay modes as discussed in the text. The two entries for D — » 7777 
and D° — > 7777' correspond to the final states that are reconstructed with two 77 — > 77 decays, or one 77 — > 77 and one 
77 — » 7r + 7r _ 7T . 



DULII LL 


qr- 

77 7T 


/ + 

77 7T 


u 

777T 


/ u 

77 7T 


7777 


1 

7777 






+ 7r - 


„' 7r + 7r - 


IVIV^ OLailSLlCS 


1 n 


1 n 

l.U 


1 n 

l.U 


1 n 

l.U 


1 n 

l.U 


1 n 

l.U 


9 n 
z . u 


9 n 
z . u 


1 n 

l.U 


9 n 


Particle Recon. & ID 


3.7 


3.9 


4.3 


4.3 


7.4/4.3 


7.5/4.5 


4.3 


4.4 


3.8 


4.0 


Ai? Selection 


2.0 


2.0 


4.0 


3.0 


3.0/2.0 


2.0/2.0 


4.0 


3.0 


3.0 


2.0 


77' Selection 




1.0 




1.0 




1.0 




1.0 




1.0 


Signal Shape 


1.0 


3.0 


1.0 


3.0 


3.0 


2.0 


2.0 


3.0 


3.0 


4.0 


Background 


0.6 


0.5 


2.7 


3.2 


1.4 


3.1 


1.6 


4.4 


2.7 


3.2 


Multiple Candidates 


1.0 


1.0 


1.0 


1.0 


1.0/3.0 


1.0/1.0 


1.3 


1.5 


1.0 


1.0 


77 — > 77 Veto 














5.0 
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FIG. 4: Sideband-subtracted invariant mass distributions for ?77r + (left) and it + tt~ (right) in the D° —± t]it + tv~ decay. 
The points are data and the histogram is a phase space (PHSP) model of the decay from Monte Carlo simulation. 

for CS decays. The topological amplitudes for these decays are given in Table fVl We first assume SE 1 = 
and using these four measured branching fractions, we fit for C", E' and cos Sc'E'- The resulting amplitudes 
and relative strong phase are shown in Table IVH and are compared to the values from CF decays, scaled by 
A = tan6»c7 = 0.2317 The values agree within about two standard deviations, and are well within the 
oft cited ^20% level of SU(3) symmetry breaking effects in D meson decays. 

We now allow for an additional singlet exchange amplitude, SE' The SU(3) representations, including 
SE', are shown in Table El Invoking SU(3) symmetry, we have T'"= XT, C = AC, E' = XE and A' = XA, 
5c"T' = $ct, Se'T' = Set, and §a'T' = Sat- The amplitudes T' , C , E' , and A' (in units of 10~ 7 GeV) are 
found to be: 
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TABLE IV: Table of branching fractions, SU(3) representations [g|, decay momenta, and invariant amplitudes for 
D — > PP Cabibbo-favored decays. The last column shows the fitted branching fractions as described in the text. 
The branching fractions, £? exp , are taken from recent CLEO measurements [lOl. Il8l. [l9l| when available, otherwise we 
use PDG values [1 lj . 



Mode 




Representation 


^cxp 


P* 


A 


Bet 








(%) 


(MeV) (10" 6 GeV) 


(%) 


D (i -> 


K-TT + 


T + E 


3.891 ±0.077 


861 


2.52 ± 0.02 


3.899 


D° 


K"ty° 




2.238 ±0.109 


860 


1.91 ±0.05 


2.208 


D° - 


K° n 


1 ^ C 


0.76 ±0.11 


772 


1.18 ±0.09 


0.76 


D° -> 


K° V ' 




1.87 ±0.28 


565 


2.16 ± 0.16 


1.95 


D+ -» 


K (, n+ 


(C + T) 


2.986 ± 0.067 


862 


1.39 ±0.02 


2.99 




K"K+ 


(C + A) 
js(T-2A) 


2.98 ±0.17 


850 


2.12 ± 0.06 


3.02 






1.58 ± 0.21 


902 


1.50 ±0.10 


1.47 








3.77 ±0.39 


743 


2.55 ±0.13 


3.61 



TABLE V: SU(3) representations [jj, measured branching fractions (from this analysis), decay momenta, and invariant 
amplitudes for selected Cabibbo-suppressed decays. The branching fraction uncertainties are obtained from the 
quadrature sum of the statistical and systematic uncertainties. 



Mode 




Representation 
(SE' = 0) 


Representation 
{SE' 0) 


(10- 4 ) 


p* A 
(MeV) (10" 7 GeV) 




r,7T [> 


^'-2E') 


^(C'-2E'-SE') 


6.4 ± 1.1 


846 3.26 ±0.28 


D° - 


Tj'ir° 




-^(C + E' + 2SE') 


8.1 ±1.6 


678 4.09 ± 0.40 


D° 


iji] 




^f(C' + SE') 


16.7 ± 1.9 


755 5.57 ±0.32 


D ^ 


Tff}' 


^ 2 (C' + 6E') 


f±{C + 6E' +7SE') 


12.6 ±2.7 


537 5.74 ± 0.61 



T' = 6.44 ±0.30, 

C = 4.73±0.39e i( - 151±2) ° 

E' = 3.89±0.28e^ 117±4 >° 

A' = 1.25±0.86e l (- 64 -")° 



= (-4.15 ± 0.38) +«(-2.25 
= (-1.76 ± 0.24) + i(3.48 ± 
= (0.55 ± 0.34) + 14 ± 



(4a) 

±0.15), (4b) 
0.29), (4c) 
0.83). (4d) 



TABLE VI: Solutions for the topological amplitudes using Cabibbo-suppressed decays. The fit results are compared 
to the values from Cabibbo-favored decays, scaled by A = tan#c = 0.2317. 



Amplitude 


Magnitude 


s CE n 




(HT 7 GeV) 






CS AxCF 


CS CF 


C 


5.8 ±0.3 4.7 ±0.4 




E 


3.5 ±0.3 3.9 ±0.3 


77 ± 7 92 ± 4 
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One may rewrite the amplitudes in Table |V] as follows [1] : 



V6A(D' 




IE' -C + SE' = (0.63 + i9.21) x 10~ 7 GeV + SE', 
-{C + E') + SE' = (-2.95 + zO.61) x 10~ 7 GeV + SE' 



(5a) 




C + SE' = (-4.15 - i2.25) x 10~ 7 GeV + SE' 



-(C + 6E') + SE' = (-2.10 + i2.66) x 10~ 7 GeV + SE' 



(5b) 



(5d) 



(5c) 



The right hand side of each of these four equations defines a vector in the complex plane, which contains 
an unknown complex offset, SE'. The left hand side of these equations defines four circles in the complex 
plane, whose radii are determined by their measured branching fractions. We thus have four constraints 
with which to solve for the real and imaginary parts of SE' . The solutions are obtained graphically from 
the common intersection point (s) of these four circles. 



Figure [5] shows Eqs. I5all5dl graphically, where we find two common intersection points, which are (the 
negative of) the allowed SE' solutions. The small solution Re(S'£") = (-0.7 ± 0.4) x 10 -7 GeV and 
Im(SE') = (-1.0 ± 0.6) x 10~ 7 GeV, is favored due to the expected OZI suppression. The larger solution, 
which is disfavored, corresponds to Re(SE') = (5.3±0.5) x 10~ 7 GeV and lm{SE') = (-3.5±0.5) x 10" 7 GeV. 
We thus find that apart from a small additional SE' contribution to CS decays, C ~ AC and E' ~ XE, and 
thus these decays respect SU(3) symmetry at the level of ~20%. We therefore find that this SU(3) flavor 
diagrammatic approach provides a reasonable description of both Cabibbo-favored and Cabibbo-suppressed 
D — ► PP decays that depend on C and E diagrams. This SU(3) topological approach may not apply equally 
well to all decays. For example, it is well known that the prediction B(D° — > tt + tt~) = B(D Q — > K + K~) is 
not realized due to SU(3) symmetry-breaking effects in the form factors and decays constants. 

In summary we report on first observations of D° — > r/'ir , D° — > rjrj, D° — > r\rj ', and D° — > r)ir + ir~ . Wc 
also find evidence for the three-body decays D + — > 777r + 7r°, D° — > r/7r + 7r~ and D + — > ij'ir+ir . We have 
analyzed D° — > 7777° , tj'tt , r\r\ and D° — > 777/ decays within the SU(3) flavor topology approach P and find 
that the color-suppressed and exchange amplitudes have magnitudes and a relative strong phase that are 
consistent with Cabibbo-favored decays. We have performed a second fit where we allow for an additional 
singlet-exchange amplitude and find a solution near zero (favored) and a larger solution, which is disfavored 
due to the OZI suppression of this process. 
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